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Abstract 

The purpose of this paper is to analyze the so-called back-off technique of the IEEE 802.11 protocol 
in broadcast mode with waiting queues. In contrast to existing models, packets arriving when a 
station (or node) is in back-off state are not discarded, but are stored in a buffer of infinite capacity. 
As in previous studies, the key point of our analysis hinges on the assumption that the time on 
the channel is viewed as a random succession of transmission slots (whose duration corresponds 
to the length of a packet) and mini-slots during which the back-off of the station is decremented. 
These events occur independently, with given probabilities. The state of a node is represented by a 
two-dimensional Markov chain in discrete-time, formed by the back-off counter and the number of 
packets at the station. Two models are proposed both of which are shown to cope reasonably well 
with the physical principles of the protocol. The stabillity (ergodicity) conditions are obtained and 
interpreted in terms of maximum throughput. Several approximations related to these models are 
also discussed. 


Keywords: Analytic model, broadcast, IEEE 802.11 protocol, Markov chain, ergodicity, 
generating function, Laplace transform, random walk, wireless LAN. 
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1 Introduction 

Several studies mil! have recently been devoted to the analysis of the IEEE 802.11 protocol, both 
in the unicast and broadcast modes. Nonetheless, to the best of our knowledge, no model has 
taken into account the possibility of having buffers to store arriving packets which, due to channel 
occupancy, can not be transmitted. When using the framework developed in min, it is impossible 
to rigorously evaluate the stability of the protocol. In this paper we use the same key assumption 
concerning the slots as in mii!, but we couple it with a Markovian analysis of nodes using the 
IEEE 802.11 back-off with an infinite buffer for the packets. We discuss how the way of sensing 
slots on the channel can give rise to two closely related models. The first one is analyzed in Section 
and considers the so-called greedy mode, while the second one in Section deals with a fair 
load situation. In both cases, ergodicity conditions and the value of the maximum throughput are 
obtained. The waiting time distribution of a packet is tackled in Section 

2 System parameters 

In contrast to the famous ALOHA protocol where the back-off does not take into account the 
activity of the channel, the back-off scheme of IEEE 802.11 monitors the channel in order to 
schedule the transmission of its pending packets. When a node receives a packet to send, it first 
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3 Model of an isolated station in a greedy mode 


selects a random back-off number, which indicates the number of mini-slots the node has to wait 
before transmitting. A node with a pending packet senses the channel, and decreases its back-off 
counter by one each time an idle mini-slot is detected. When its back-off counter reaches the value 
zero, the node transmits the packet. 

The following parameters are ubiquitous in the analytic studies of IEEE 802.11 broadcast: 

• (T, the duration of a mini-slot (MS), which is the time needed for a station to sense whether 
the channel is busy. 

• The backoff window size W, expressed as a multiple of a. 

• T, the packet duration, a priori larger than a. 

The three above parameters are really essential in the IEEE 802.11 back-off scheme. They lead to 
the main assumption that is usually made in the analysis of this scheme (see [2 H]), stating that 
the channel consists in a random succession of full slots (dedicated to packet transmission) and 
mini-slots (short slots during which the channel is idle). To complete the description of the model, 
we also need to introduce the following three additional quantities. 

• A, the Poisson arrival rate at a station. 

• M, the number of stations in the network. 

• r, the probability that a station will transmit a packet, stressing that this can take place 
only when its backoff counter has reached zero. We will return to this parameter later. 


3 Model of an isolated station in a greedy mode 


This section is devoted to the analysis of the basic component of the network, namely a single 
station operating in compliance with some of the principles of the IEEE 801.11 protocol. 


The channel is supposed to be sampled at discrete time instants Zi, so that it is sensed during 
consecutive time intervals Zi — Zi_i,i > 1, which are equal either to T (a normal slot) or to a 
(a mini-slot (see Section 3.11. In some sense, we have thus introduced two embedded time-scales. 


In general the selection between these two types of slots is random: a slot is a normal slot with 
probability r, and a mini-slot with probability 1 — r. Nonetheless the node is said to be greedy, 
since, when the back-off reaches 0, a transmission takes place in a slot which is then necessarily a 
normal slot of length T. In particular, this implies that a sequence of slots is not exactly obtained 
as a pure outcome of repeated Bernoulli trials. 


Now, starting from the basic one-dimensional models proposed in mm , we construct a two- 
dimensional Markov chain, taking into account the possibility of waiting queues at the stations. 


An arbitrary station (or node) will be represented by a stochastic process, which is a two- 
dimensional random walk {Kz^, Nz-)i>i, where Zi is an embedded increasing sequence of discrete 
random times at which the chain is observed. Indeed, at time Zi, Kz^ stands for the value of the 
backoff counter and Nzi is the number of packets in the waiting queue at a station. The Markovian 
evolution of {Kzi,Nzi) is described in detail in the next section. 


3.1 Dynamics 

Taking into account the definitions given in Sectionthe sequence Zi,i > \ va the so-called greedy 
mode satisfies the recursive stochastic relationship 

Zi+l = Zi -\- Ai[l[Kz->o} + ^{Kzi=Nzi=0}] + ^l{ifzi=0.Arzi>0}i (3-1) 

with 

Ai = TI1 {bj-I- cr(l — I1{B;}), (3.2) 

where Il{ j denotes the indicator function and Bi is the event {channel busy} at time Zi. From the 
assumptions made above, {Il{Bi})* > 1} form a sequence of independent identically distributed 
random varaibles with V{Bi) = r. 
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Without entering a quite standard formalism, we just state that the underlying probability space 
for the random walk {Kz^, Nzi)i>i is obtained by combining the external Poisson arrival process 
and the sequence of independent Bernoulli trials related to the Si’s. 


Then we are in a position to write Kolmogorov’s forward equations for the Markov chain , Nz^). 
The reader will see that the greediness appearing in the title of the section stems from the last 
term in the right-hand side member of (3.1). 


Let us define the conditional probability p{k,n\Zi) = P{Kzi = k,Nzi = n;Zi). Then, for all 
k > 0,n > 1, 


p{k,n]Z^+T) = ll{fc>o}ll{B,} ^^— -p{k,n- j;Z,) 


■E 


W + l^ J 

J=o ■’ 


j=o 

-j - p{0, n + l-j,zy 


1 e-^'^(AT)" 


U 




W+1 nl 
p{k,n;Zi + a) = ll{fc<w}(l - l{Bd)E 


p{0,0;Zy 

(Aay 




J=o 


J- 


p(k + l,n - j; Zy 


1 e“^‘^fAcr')” 


(3.3) 


(3.4) 


These equations deserve some explanation. 


The first term of the right-hand member of (3.3) corresponds to the transmission of a packet 


by another station when the current station is in back-off, noting that arrivals may occur 
during this slot. The second term corresponds to a packet transmission by the current station 
and the choice of a back-off number for the next pending packet; new arrivals may still occur. 
The last term concerns arrivals at an empty current station during the transmission of a 
packet by the network; a random back-off is selected for the first witing packet (if any). 


It is worth remarking that for k = 0 the first term in the right-hand member of (3.3) 


does not exist, since in this case it is impossible to reach the state (0, n, Zy from any state 


{0, j, Zi_y, j < n. Also, the second term in (3.3) indicates that a transmission takes place 
when the backoff counter equals 0. 


In (3.4), the first term of the right-hand member corresponds to the current station 


decrementing its back-off counter; new arrivals may occur during the mini-slot, while the 
second term takes into account the new arrival during a mini-slot in the current station 
which is idle; if at least one packet has arrived a random back-off number is selected for this 
packet. 


On the other hand, the following boundary equations hold for k = 0 and for all Ziy > 1. 


(p{k,o;zy = o, yk>i, 

I p(0, 0; Z, + T)= e-^^ [p(0,1; Zy + p(0,0; Zy] , (3.5) 

[p(0, 0; Zi + a) = e“^‘^(l - 11 {b,}) p(0, 0; Zy, 


where p(0,0, Zy is the empty (or idle) state, through which the system exhibits some regeneration 
properties. 


In the following, we focus on the steady state behaviour of system ( |3.3[ )-(3.4)-(3.5), which is simply 
obtained by letting i — >■ oo, since in this case —>■ oo almost surely (and even surely!). 


Setting 


p{k,n) = lim E[p(fc, n; Z^)], 

i—>-oo 
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3 Model of an isolated station in a greedy mode 


we obtain, for all n > 1, 0 < fc < W, 


and 


p(k,n) = - J) 

j=o 

. ^e-^^ixry ,, ,, 

+ l{fc>o} r 2 _^ -^- p[k, n- j) 


j=o 

” ^-AT/ 


+ 


+ 


1 {XTy 


w- 




f=o 


J'- 


p(0,n + l-j) 


1 


W + ll 


{1 — r)e re ^^(AT) 


n! 


p(0,0), 


I p{k, 0) = 0, VI < fc < W, 

e~^^p{0, 1) = p{0, 0) [l — re~^^ ~ (1 ~ r)e~^'^]. 


(3.6) 


(3.7) 


Let us observe that here r is an exogenous parameter, but later in Sections [3.4[ |4.1| it will depend 
on the state of the system. 


3.2 Formal solution of the steady-state equations 

Let us introduce the generating functions 

Fk{x) = ^p{k,n)x'^, V0<fc<LL, 

n>0 


which will be sought to be analytic in the unit disc T) = {x : \x\ < 1} and continuous on the 
boundary \x\ = 1. 

Lemma 3.1. The generating functions Fk{x),k > 0, satisfy the linear system 


( a{x)Fi{x) = Fq{x) + C{x), 

< a{x)Fk+i{x) = b{x)Fk{x) + C{x), VI < /c < VL — 1, 
\^b{x)Fwix) + C{x) = 0, 


where 


and 


C{x) 


r?b(a^) = i^o(a;)-p(0,0) 
< a{x) = (1 - 
[b{x) = 


_g-Ar(l-a:) 
(W + l)x 


Fo{x) + 


p{0,0) 
W + 1 




(1 - 


(3.8) 


(3.9) 


(3.10) 


Proof. System (3.8) is easily obtained after mutiplying (3.61 by a;", summing from n > 1 onwards, 
and taking into account the first boundary condition in (3.7). The final form of C{x) in ( |3.10 ) is 
derived from the second boundary condition in (3.7). More detail is given in Appendix [A| ■ 


Proposition 3.2. Letting 


. def b{x) 
u{x) = 

a{x) 


the functions Ff^^x), 1 < k < W, are given by the relations 


Fk{x) 


Fo{x) ^{x) — u^{x) 
a(x) _ 1 — u^'^^{x) 


(3.11) 
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and Fg (x) has the explicit form 


Fo(x)=p(0,0) 


Q{x) 

R{x) ’ 


where 


R{x) = 


e {W + l)u^{x){l — u{x)) 


Q{x) = 




1 - u^+i(a;) 


— T-— (1 — 


(W + l)u^(a:)(l — u{x)) 


(3.12) 

(3.13) 

(3.14) 


l-u^+i(x) 

Proof. The result is obtained in two steps. First a straightforward calculus exploiting 
system (3.8) yields (3.11[). Then, instantiating fc = IF in equation (3.11) and comparing it with 


the third relation in (|3.8[), we get 


_- 1 _ 

u^ix)Fo{x) + C{x) ■ ^ ^ f = 0, 
1 — u[x) 


which, by using (3.9), leads to the final equations (3.12),(3.13),(3.14). ■ 

The only remaining unknown p(0, 0) will now be obtained from the normalization condition 

X! P{k,n) = l. 

k>0,n>0 

Proposition 3.3. Under ergodicity conditions, the probability p(0,0) that the system is idle is given 
by the formula 

p(0,0)= 


1 _ \4 _L >^W(B-A) 
i AA+ 2(l-r) 


where 


'A =(l-r)(r-a)+ ^t-^+y-] . 


p _ 'T' \ VK[rT+(l —r)cr] 

^ ^ 2(l-r) 


(3.15) 


(3.16) 


Proof. By ( |3.11[ ), since m( 1) = 1 and a(l) = 1 — r, we get 

w 

w - 

ft—i 

that is 


Fin I + ^ 1 


(-) 

Hence, instantiating x = 1 in equations ( 3.12| ),(3.13)-(3.14), using the relation 

Q{x) Q'(i) 

x™F(x) F'(l)’ 

since Q(l) = F(l) = 0, where Q'(l) and F'(l) are the respective derivatives of Q{x) and R{x) at 
X = 1, from equation (3.17) we obtain 


p(0,0) 

Setting for the present calculation 

fix) = 


1 + 


W 


] Q'W 


w 


= 1 . 


2(l-r)yF'(l) 2(1-r)_ 

(fF + l)it^(x)(l-u(x)) (fF + l)u'^(x) 


1 - u'^+i(x) 


U*(x) 


0<i<W 
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3 Model of an isolated station in a greedy mode 


we have 


/'(I) = 


Hence, by (3.13) and (3.14) 


R'{1) = _l + Ar-/'(l) = -1 + AT + 


Wu'{l) AfH[rr + (1 - r)cr] 
2 “ 2(1-r) ' 

\W[rT + (1 — r)a] 


2(1-r) 


g'(l) = -1 + (1 - r)A(T -a)- /'(I) = -! + (!- r)A(T - a) + 


\W[rT + (1 — r)a] 
2 ( 1^0 ’ 


which by an elementary computation yields equation (3.15). 


Remark 3.4. As we shall see in Theorem |3.5[ the system is ergodic if and only if 

A'(l) < 0, 

in which case the following inequalities hold 

'Q'{1) < R'{1) < 0, 

Q'(l) 


^'(1) 


> 1 . 


Two limit cases can also be checked. 

• A —?► 0 ^p(0,0)—>■!. 

• It should be observed that Q'{1) = R'{^) if and only if r = cr = 0, which corresponds to the 
trivial situation p(0,0) = 1. 

□ 

We bear in mind that functions Fk{x),k > 0, are sought to be analytic in the unit disc. Hence 
ergodicity conditions will be obtained from (3.12) by studying the possible zeros of R{x) in the 
closed unit disk V, noting en passant that Q{x) and R{x) have no common root in V, but at x = 1. 
This is the objective of the next theorem. 


3.3 Stability condition of the greedy model 


With the notation of Proposition |3.3[ we have the following 

Theorem 3.5. The system is ergodic if and only if R'{!) = XB — 1 < 0, that is 


A < 


1 


T 


1 + 


rW 

2(l-r) 


Wa 

2 


(3.18) 


In other words, the ergodicity region in the parameter space {A,(t, T, r} is delimited by the surface 
S{X,a,T) expressed by the equation 


Proof. Let 


\[{Wr + 2(1 - r))T + W{1 - r)a] = 2(1 - r). 


-AT(l-a;)ri V 

g{x) = - - , ^ ^ hix) = {W + l)xu'^{x), 


1 — u(x) 

one can see at once that, for all x yf 1, the zeros of i?(x) in V coincide with the roots of 

g{x) - h{x) = 0. 


From (3.9) 


i{x) = 


gA(T(l-a:) _ ^g-AT(l-a;)j 


(1-r) 
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so that, \/x G V, 


\ — I 

|u(a:)| = -^|1 - > 1, 


1 — r 


which yields 


1 — 

1 — u(x) 


|1 + u{x) H-h U^{x)\ < {W + l)\u^{x)\. 


Therefore, on the unit circle \x\ = l,a: 7 ^ 1, we have |g(a;)| < |/i(a;)|- Since the functions g{x) and 
h{x) are analytic in I?, we will be in a position to apply the well-known Rouche’s theorem (see e.g. 
0 ), after carrying out an adequate modification of the unit circle C around the point x = \. Since 


g{x) - h{x) = -^--- R{x), 

1 — u[x) 


we have 

g{x) - h{x) = {W + - 1) -b o{{x - 1)^). 

Let us assume R'{!) < 0. We construct a contour Cg consisting essentially of the unit circle C 
continuously distorted around a: = 1 by a small notch (an arc of a circle of radius e), keeping the 
point 1 inside the domain bounded by Then, choosing e sufficiently small, we have 


\g{x)\ < \h{x)\, VxGCe. 


Consequently, by Rouche’s theorem, functions h{x) and h{x) — g(x) have the same number of zeros 
inside V^, i.e. exactly one located at x = 1, since h{x) vanishes only at x = 0. We conclude that 
Fo(x) given by (3.12) is analytic in V and the system is ergodic. 

Conversely, if R'{1) > 0, we can construct a contour that does include the point x = 1, and 
the same argument shows that R(x) has a root different from 1 inside the unit disk I?, and in this 
case Fo(x) is not analytic in V. ■ 


We observe that, on the one hand, for r tending towards 1, A tends towards 0. The activity 
outside the current station is prominent, which means that no bandwidth is left for the current 
station. On the other hand, if r —)■ 0, then A can reach the value Amax = T+Waj’i ■ 
the current station can consume all the bandwidth, the term being the overhead due to the 
back-off scheme. 


3.4 Dynamics of a network of greedy broadcasting stations 

Consider a network consisting of M identical stations operating in a broadcast mode together 
with a tagged station. The key feature of the protocol IEEE 802.11 —namely a transmission is 
permitted only when the backoff counter is at zero— can be reflected by the simple equality 

r = ^p(0,n). (3.19) 

n>l 

Hence, from the point of view of the tagged station, the M stations represent the outside world 
and the event { Channel husy\ takes place with probability 

r = 1 - (1 - M > 0. (3.20) 


The problem is now to recast the ergodicity condition ( |3.18| ) given in Theorem 3.5 in the light of 
equations (3.19), (3.20). According to Propositionand Remark [3)4{ r defined by (3.19) satisfies 


r = Fo(l)-p( 0 , 0 )=p( 0 , 0 ) 


Q'W 

R'{1) 


-1 =p(0,0) 


\{B-A) 
l-XB '■ 


that is, by ( 3.15| ) and (3.16), 
A(B - A) 

. ... XW{B-A) 

i A^l-I- 2(l-r) 


A[rT -b (1 — r)cr] 


A[rT -b (1 — r)cr] 


1 _ \4 _L >'W{,B-A) 
1 2{l-r) 


1 - XT + X[rT + {1 - r)ay 


T = 


(3.21) 




































3 Model of an isolated station in a greedy mode 


Observe at once that the denominator in the last equality of (3.211 necessarily imposes 

AT < 1, 


(3.22) 


with an obvious interpretation. Henceforth condition (3.22) will be assumed in the rest of this 
section. 


So, (3.20) and (3.21) yield the following fixed point equation in r, written as 

1-AT 


(l-r)i/^ = 


On the other hand, setting 


1 - A(1 - r){T-a)' 


z= (l-r)i/“, 


one can also view (3.23) as a polynomial equation in z, namely 

P{z) = A(T - a)z^+^ - z + (1 - AT) = 0. 


(3.23) 


(3.24) 


(3.25) 


Lemma 3.6. Under condition (3.22), equation (3.25) viewed as an equation in z has only one root 
in the the real interval [0,1]. 

Proof. Under the assumption AT < 1, we have 


P(0)>0, T(1)<0, 


and two cases must be considered. 


(i) T < a (purely academic situation). Then, on the interval [0,1], the derivative P'{z) in (3.24) 
is negative, and P{z) vanishes exactly once. 

(ii) T > a (the real world!). On the unit circle |z| = 1, we have 

|z + l-ATI > ||z| - (1-AT)| > AT> \X{T - a)z^+^\ = A(T-ct), 

and the result follows by an immediate application of Rouche’s theorem. 

The proof of the lemma is concluded. ■ 

We adopt the following figures A = 0.05, T = 1, tr = 0.05, W = 31, and we vary M. We compute 
r in our model and we compare it with the value of r obtained in the model without any buffer 
of [3]. In Figure 3.1 we observe that r is larger in our model compared to the model without 
any buffer for the same input load. This is because the model without buffer drops packets when 
one packet arrives during the back-off period of a preceding packet. In Figure 3.2 we compare 


the throughput without collision of our model with the throughput without collision of the model 
without buffer. For the same reason, the throughput without collision of our model is higher. 

We are now in a position to compute the maximum throughput of the system under the stability 
condition (3.18). This is the subject of the following and final theorem. 

Theorem 3.7. 

1 The greedy broadcast network is ergodic if and only if 

2z^+i > lU(l-z), 


(3.26) 


where z defined by (3.24) is the unique real root of equation (3.25) in the interval [0,1]. 

2 When the system is ergodic, the maximum achievable throughput of the system takes one of 
the two equivalent forms 


Xmax — 


I — U 


I — U 


T(l-u^^+l)-hcru^+l 
where u is the unique root on [0,1] of the equation 


W(a-T)(l-u) 

2 


T- 


(3.27) 


= IU(l-u). 


(3.28) 
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Fig. 3.1: Computation of r with and without any buffer, letting vary the number of 
stations. 


Proof. Rewrite B given by (3.161 as 


B = T- 


W[{l-r)[<7 -T) + T] 
2 (1-r) 


which, using (3.251, yields the pleasant factorized form 

1 - AS = (1 - XT) 


. Wjl-z) 

22M+1 


Hence, the necessary and sufficient ergodicity condition XB — 1 < 0 comes down to 

1- '"'(!-») >0 
2s'W+> ’ 


which is exactly (3.261. Finally, the quantity Amax is simply obtained by saturating inequality 


(3.26), which gives (3.28) and then (3.27). 


We exploit formulas ( |3.27[ ) and ( |3.28[ ). We still have T = 1, cr = 0.05, W = 31, but A = 0.05, 
and we vary M between 1 and 100. We plot the maximum throughput of one station and the 
cumulated throughput of all the stations which we call the network offered load. This is shown 
in Figure [33| where we vary the number M of stations in the network. One can observe that the 
network offered load can be larger than 1. This can be explained by the fact that two stations (or 
more) can send a packet during the same slot. The network offered load is the maximum total load 
the network can submit while remaining stable. For a given number f stations, it is also possible 
to study the maximum network throughput (without collision) when we vary W. This study is 
presented in Figure [3l4| 


4 A fair load model 

Here a station does not impose an extra load, as this was the case in the greedy model, and the 
sequence of times Zi,i > 0 is given by recursive scheme 
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4 A fair load model 



Fig. 3.2 : Network throughput without collision versus the number of stations. 


where is still given by (3.21. A station with a pending packet and a back-off counter at 0 may 


send its packet with probability r or return to back-off with probability 1 — r. Thus, in contrast 
to Section slots on the channel are selected at random and form a renewal process: a full slot 
(length T) with probability r, or a mini-slot (length a) with probability 1 —r, independently of the 
state of the station. So, r is considered as an exogenous parameter. But if the station is able to 
sense the channel, one could think of r being a function of the state, possibly estimated via some 
adaptive scheme like for example in [5] for Aloha-type systems. 

Mutatis mutandis with respect to Section |3.1[ one sees that the stationary distribution of a fair 
loaded station satisfies the following Kolmogorov’s equations, for all n >, 0 < fc < FF, 


q(k,n) = a{fe<vy}(l-r)^ 


o —Act 


(Aa)* 


q{k + l,n-i) + ll{fc>o} r ^ 


,-AT 


{xry 


q{k, n — i) 


+ 


1 


W + 1 


•E 


2 = 0 

e-^'^{XTy 




z: 


9(0,n-bl-f)-f (l-r)^ 


2 = 0 
A(T 


(Act)’ 




2=0 


2 ! 


q{0,n-i) 


re“^^fATl" 

+ 9 ( 0 , 0 ), 


{W + l)n\ 


(4.1) 


and 


q{k,0) = 0, VI < fc < W, 

re~^^q{0,1 ) = q{0, 0) [l — ~ (1 ~ r)e“^‘^]. 


(4.2) 


The first line of equation (4.1) corresponds to the transition when the back-off counter is greater 
than or equal to 1. If the current slot is a mini-slot (with probability 1 — r) the back-off counter is 
decremented, otherwise if it is a normal slot (with probability r) and the back-off counter remains 
unchanged. In both cases new arrivals may occur and increment the number of pending packets. 
The second line of equation (4.1) deals with the transition when the back-off counter is at 0. With 
probability r, there is a transmission with a new draw of the back-off (if any pending packets) and 
new arrivals may occur. With probability 1 —r the current slot is a mini-slot and the packet returns 
to back-off, and there still may be new arrivals. The third line of equation (4.1) corresponds to 
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Fig. 3.3 : Network throughput versus the number of stations M. 


the transition from the idle state of the station to a state where a back-off value is selected with 
probability 1/{W + \). 

Let us introduce as before the generating functions 

Gkix) = ''^q{k,n)x'^, 0<k<W, 

n>0 


sought to be analytic in the unit disc V = {x : \x\ < 1} and continuous on the boundary \x\ = 1. 
Then the following results hold, the explicit proofs of which are omitted, as they mimic the 
arguments of Section |T2| Some hints are given in Appendix |B| 

Lemma 4.1. The generating functions Gk{x),k > 0, satisfy the linear system 


{ a{x)Gi{x) = Go{x) + D{x), 

a{x)Gk+i{x) = b{x)Gk{x) + D{x), VI < fc < Fh — 1, (4.3) 

b{x)Gw{x) + D{x) = 0, 


where 


a{x),b{x) being given by (3.9|, and 
-Go{x) 


D{x) = 


re 


-AT(l-x) 


{W + 1) 

Proposition 4.2. Letting 


Go(x) = Go{x) - g(0,0). 


-F (1 - 


] ^ g(o,o) 


J W + 1 

\x ) 


u{.) a 


(x)’ 


the functions Gk{x), 1 <k < W, are given by the relations 


Gk{x) 


Go{x) ^{x) — u^{x) 
a{x) _ 1 — u^'^^{x) 


(4.4) 
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4 A fair load model 
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Fig. 3.4: Optimized network throughput (without collision) and corresponding value of 
W versus the number of stations M. 


and Gq (x) has the explicit form 


Go{x) = <7(0,0) 


Qjx) 
R{x) ’ 


where 


R{x) = 
Q{x) = 


.xaii-y) _ {W + l)u'^{x){l-uix)) 
^ ’ l-u^+^{x) 


1 + rf- - ] 


\x 

/ 


(W + l)u^(x)(l — u{x)) 


l — v}^+^{x) 

Proof. Along arguments quite similar to those used in Section |3.2[ one can show 


...^ 1 _ 

u'^{x)Go{x) + D{x) —---— = 0, 

1 — u{x) 


which directly yields (4.5), (4.6|, (4.7). 

Theorem 4.3. The fair load station is ergodic if and only if R (1) < 0, that is 

r(l - r) 


A < 


In this case. 


9 ( 0 , 0 ) = 1 - 


[1 - r + 14^/2] [rT + (1 - r)a]' 
X[rT + (1 — r)a] 


1 + 


w 


2(l-r) 


(4.5) 

(4.6) 

(4.7) 


(4.8) 


(4.9) 


Proof. As in Proposition 3.3 by using (4.4) and (4.5), the normalization condition becomes 



























































4.1 Dynamics of a network of fair broadcasting stations 
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whence, after some algebra 


g(0,0) 


1 


w 

2(l-r) 


Q (1) 
.it'(l) 



+ 


Q (1) 

«'(i) 


-i?'(l) 


(4.10) 


which is exactly equivalent to (4.9). One can remark that when r 0, the system remains ergodic 
if and only if 


A 

r 


1 

a(l + W/2)’ 


in which case g(0, 0) \ 0. 


4.1 Dynamics of a network of fair broadcasting stations 

We proceed as in Section |3^ setting 


with 

r = ^(?(0,n). 

n>l 

Then the following result similar to Theorem |3 . 7| holds. 

Theorem 4.4. The fair broadcast network is ergodic if and only if 

A Xrnax'! 


where the maximum achievable throughput \max takes the form 

A - 

/'max — 


W(7{l — u) 

u(2+W)-W 


T- 


where u is the unique root in [0,1] of (3.28), that is 

= W{l-u). 


(4.11) 

(4.12) 


(4.13) 


Proof. 

We outline the arguments, as they mimic those of Theorem |3.7| The algebra is even simpler. 
Indeed, by (4.12) and (4.10) one can write 

T=,(0.0) (gai _ W g. d) - , x[rT _ 

\i?(l) / ^ ^ 


so that, using (4.8) 


lim r = 


1 


A—f An 


1 + 


W 

2(l-r) 


Then, putting 1 — r = , it appears that u satisfies exactly (3.28), and we obtain 

Amax — ~ 


1 


1 — U 


1 + 


w 

2(l-r 


T- 


( 1 - 0 ° 


T + 


P-'r)o 


which with (4.11) and (4.14) yield (4.13). 


(4.14) 


It is clearly worth comparing the two maxima Amax and Amax given respectively by (3.27) and 
(4.13). To this end, it suffices to check that the quantity 

a — T a 


6 = 


1 


u{2 + W)-W 2 


a-T- 


u(l — u) 


M 
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5 About the waiting time distribution in the greedy model 


is clearly negative as m < 1. Hence 


A„ 


< A„ 


We exploit the formula giving the maximum station throughput. We still have A = 0.01, T = 1, 
(T = 0.05, W = 31 and we vary M between 1 and 100. We plot the maximum throughput of one 
station for the fair and greedy models. This is shown in Figure |4^ where we vary the number M 
of stations in the network. As foreseen, the maximum achievable throughput is lower in the fair 
model compared to the greedy model. However, in both models u is given by the unique equation 


(3.28), so that the station transmission rates coincide, as well as the network throughputs without 


collision. 



Fig. 4.1 : Maximum station throughput versus the number of stations M in the fair and 
greedy models. 


5 About the waiting time distribution in the greedy model 


From now on, the system will be assumed to evolve at steady state. We intend to compute the 
stationary distribution of the virtual waiting time of a packet arriving at an arbitrary epoch Zi. 
For an arbitrary probability distribution F{.), denote by F* its Laplace-Stieltjes transform (LST) 


pOO 

F*{s) =* / e~’^*dF{t)dt, 

Jo 


3?(s) > 0. 


Setting 


/(s) 


(1 - 

1 — ’ 


n(s) 


e-*^[/'^+i(s)- 1] 

(W + l)(/(s)-l)’ 


(5.1) 


we have the following result. 

Theorem 5.1. The Laplace transform of the stationary virtual waiting time distribution of a packet 
arriving at time Zi has the form 


w 

V'*(s) = J2f{s)Fkiv{s)), 




(5.2) 


In particular, the expectation can he computed from formulas (3.11), (3.12). 
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Proof. 


Let A(k, n; Zf) be the conditional waiting time of a packet which at time Zi sees the system in the 
state {Kzi.Nzi) = {k,n). With the notation of Section 3.1 we can write 


A{k,n;Z,_i) = (1 - lL{B._^})[cr + A(fc - l,n;Zi)] + 1L{b._j}[T + A{k,n]Zi)], V/c > 1, (5.3) 

remembering the transmission discipline at the buffer is First-In-First-Out (FIFO), without 
forgetting the boundary equation for k = 0, namely 

w 

A{0,n-,Z,_i) = T + Y,Mu=(}Aie,n;Z,)], Vk > 1, (5.4) 

£=0 


where U is random variable uniformly distributed on the integers 0,1,..., IF. 

Let tp*{s] k,n) be the LST corresponding to the distribution of A{k,n; Zi). Then (5.31 implies 

Lp*_i{s] k,n) = (1 — r)e~^'^ip*{s; k — l,n) + re~^'^ip* (s; k,n), Vfc > 1, 


whence, letting * —>■ oo and tp* = limi_,,oo p* and using the definition (5.11, 

(p*{s; k, n) = f^{s)^*{s; 0, n). 


(5.5) 


Similarly, from (5.4), the LST (p*(s;0,n) satisfies the following stationary recursive relationship 
holds. 


(p*{s;0,n) = 


,-sT 


IF + 1 


w 


1 ) = 


e=o 


r+Hs)-i 


e ®^(/7*(s; 0, n — 1) 
IF +1 ■ 


(5.6) 


Setting (/j*(s; 0,0) = 1, we get from (5.11, (5.5), (5.6), 

ip*{s-,k,n) = f^{s)v{s)'' 

and the sought waiting time distribution is given by 

W oo 


= '^^Pik,n)ip*{s]k,n) 


k—0n—0 


which reduces immediately to (5.2). 


Remark 5.2. To compute the waiting-time distribution of a packet at arrival epoch, one has to 
switch from a discrete-time model to a continous-time one. We leave this question as an exercise, 
just quoting that it can be solved in two steps. 


1 First, find the distribution of the stationary Markov chain {Kt,Nt) from {Kzi, Nz^), for 
instance by applying a Palm probability inversion formula to the marked point process Zi 
with marks [Kz^, NzJ). But a rather direct (although somehow tedious) is to proceed 
as in the derivation of residual life times for renewal processes: fix the unique interval 
[Zi(T-), containing the arrival instant r, which allows to derive the number of new 

packets arrived between and + r. Here the stationary sequence Zi satisfies 

lim E(Zi_|_i — Zi) = [r + a(l — r)] T + (1 — r)(l — a) a, 


where a = Fo(l) — p(0, 0). 

2 Then proceed as in the proof of Theorem |5.1| 


□ 
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B Equation (4.3| of the fair load model model 


6 Partial conclusion 

We have proposed two models for a CSMA protocol using the back-off of IEEE 802.11, one being 
very close to the real protocol for broadcast packets. In contrast to existing works, we assumed 
that nodes can store arriving packets. With a classical formalism of normal slot and mini-slots, 
it was possible to compute several important characteristic parameters: station transmission rate, 
stability conditions, delay for a random packet. 

It is worth to emphasize that these protocols can stabilize without external control, in contrast to 
many famous pioneers like ALOHA, CSMA, etc. In this respect, the window-backoff mechanism 
exhibits some common features with the so-called tree (or stack) algorithms, which are stable 
provided that the external load remains below some critical value. 

The numerical results obtained so far in this paper seem to match the performance of the lEEEE 
802.11 broadcast protocol. More simulations (left as future work) should very certainly confirm 
this matching. At last, studying the IEEE 802.11 protocol for point-to-point traffic under the 
same hypothesis of infinite buffer, via the same mathematical approach, might be an interesting 
follow-up of this paper. 
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APPENDIX 


A Equation ( 3 . 8 ) of the greedy model 


Next is the intermediate step leading to system (3.8). Indeed, mutiplying (3.6) by a;", summing 


from n > 1 onwards, and taking into account the first boundary condition in (3.7), namely 
p{k,0) = 0,Vfc > 1, we get 

Fk{x) = ll{fc<w}(l - r)e~^'^^'^~^'>Fk+i{x) -f Fk{x) 


1 


(W + l)x 

p(0,0) 


e-^^(i""^)(Fo(a:) -p(0,0)) -p(0, l)e-^^a; 

(g-AT(l-:r) _ g-AT) _ ^-Aa^ 


(W+l) 
which yields in particular 

a{x)Fk+i{x) = b{y)Fk{x) F C(x), VI < k < W - 1 


whence (3.8) follows easily. 


B Equation ( 4 . 3 ) of the fair load model model 


from (^4.1|,(4.2), 


Just as in the case of the greedy model, (4.3) is obtained by means of the following system, derived 


Gfc(a:) — Il{fc=o}'7(Oi0) + l{fe<iv}(l ~ ^^Gk+i(x) + ^^Gk(x) 

e-A^(i-")(Go(cr) - q(0,0)) - p(0,l)e-^^x 


{W Fl)x L 
I — r 


{W + l) L 


^-^<r(i-+Go{x) - q{{),())e 


— \a 


+ 


r( 7 ( 0 , 0)(e — e 

W + l 
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